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CHARLES  R.  COOKE 
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ABSTRACT 


This  report  describes  the  calibration  and  an  error  analysis  of  the 
6-component  50,000-pound  Gilmore  Rocket  Motor  Thrust  Stand  lo.  sted 
on  Pad  2,  Air  Force  Rocket  Propulsion  Laboratory,  Solid  Test  Area 
(1-32).  The  purpose  of  this  report  is  to  establish  stand  accuracy  and 
reproducibility  for  a  calibrated  static  test  loading. 

Errors  considered  primarily  resulted  from  thrust  stand  misalignment, 
propellant  mass  flow  of  the  motor,  and  cross-axis  interactions.  Dynamic 
effects  resulting  from  vibration  were  not  included  in  the  detailed  analysis 
but  are  briefly  discussed. 
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SECTION  I 


INTRODUCTION 

A  six-component  thrust  stand  was  designed  by  Gilmore  Induatriea  for 
me  with  an  84-inch*d>amc,ter  teat  motor.  The  B4-inch-diamcter  burn 
surface  and  typical  uncured  propellant  formulationa  produce  thruat  levels 
up  to  50,  000  pounda,  and  typical  maximum  burn  duratlona  of  60  aeconda. 
The  weight  of  a  fully  loaded  motor  ia  approximately  50,000  pounds.  With 
the  thruat  stand  framework,  the  maximum  initial  weight  ia  approximately 
60,  000  pounda. 


¥ 


The  Gilmore  thruat  atand  helps  fulfill  thruat  vector  control  (TVC) 
ayetem  teating  requirements.  The  importance  of  accurately  determining 
the  thruat  vector  produced  by  a  TVC  system  led  to  the  error  nalysia 
that  was  performed  on  the  Gilmore  thruat  stand.  The  error  .nalyeia 
serves  to  estimate  the  expected  accuracy  of  the  measured  thrust  and  aide 
force  vectors  produced  by  a  TVC  a>atem. 


Conventional  assembly  and  calibration  techniques  were  employed  on 
the  atand.  In  addition,  a  simulated  firing  loading  condition  (subsequently 
referred  to  as  the  "verification"  teat)  of  the  atand  waa  created  to  evaluate 
the  accuracy  of  the  atand.  A  study  of  this  data  waa  conducted  to  determine 
an  estimate  of  uncertainty  to  assign  to  measurements  made  with  this 
thrust  stand.  This  report  contains  the  methods,  assumptions,  and  results 
of  that  study. 
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SECTION  II 


THRUST  STAND  EVALUATION 

A.  THRUST  STAND  DESCRIPTION 

The  Gilmore  stand  for  the  84-inch-diameter  Char  motor  is  a  vertical 
six-component  thrust  stand  capable  of  measuring  50,  000  pounds  of  thrust 
and  side  forces  up  to  10,  000  pounds.  The  stand  has  an  axial  load  cell  and 
five  side  force  load  cells.  Three  10,  000-pound  load  cells  are  connected  to 
the  aft  mounting  ring.  The  remaining  two  10,  000-pound  load  cells  are 
connected  to  the  bottom  of  the  thrust  spider  next  to  the  axial  load  cell. 

The  load  cell  designated  with  the  subscript  C  (Figure  1)  are  part  of  the 
calibration  system  and  are  only  employed  during  calibration. 

With  the  geometry  illustrated  in  Figure  1,  a  Y  side  force  of  5,  000 
pounds  applied  at  the  top  will  theoretically  produce  a  force  of  47  96  pounds 
in  the  upper  load  cells  and  204  pounds  in  the  lower  load  cells.  (This 
example  illustrates  the  inadequacy  of  load  cell  error  estimates  based  on 
total  measurement  range).  The  four  calibration  cells  are  used  one  at  a 
time  to  apply  known  loads  near  the  individual  load  cells  and  develop  a 
correlation  between  cell  output  and  applied  loads  as  shown  in  Figure  2.  In 
the  verification  test,  load  distributions  similar  to  firing  conditions  are 
achieved  by  introducing  loads  into  the  thrust  stand  motor  assembly  with 
a  secondary  standard  load  cell  installed  into  the  fixture  as  shown  in 
Figure  3.  The  verification  effort  permits  the  correlation  of  load  readings 
of  the  individual  lead  cells  with  induced  loads  in  all  the  side  load  cells 
(axial  force  was  applied  separately). 

The  ideal  relationship  between  the  thrust  forces  produced  during  the 
rocket  firing  and  the  individual  load  cell  readings  can  be  determined 
from  a  forces -and-moments  balance  analysis.  Appendix  I  contains  this 
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analysis  and  the  results  are  repeated  here  to  illustrate  the  operation  of  the 
test  stand.  The  geometry  of  the  load  cells  are  shown  in  Figure  1.  Using 
the  nomenclature  of  Figure  1  the  side  forces  Fx  and  F  are  given  by: 


'Fxl"Fx3'Fx2 


+  F 


y»  y2 


The  Z  axis  thrust  magnitude  F  is  given  directly  by  the  Z  load  cell.  The 

2  /2  2  2\l/2 

total  thrust  force  vector  (T  )  is:  T  =  (F  +  F  +  F  ) 

The  thrust  vector  Ty  is,  in  general,  at  an  angle  $  from  the  motor  center- 
line  and  an  angle  $  from  the  X  axis  in  the  X-Y  plane,  measured  clockwise. 
These  angles  are  defined  by: 


$  =  arc  tan  Fz/ Ty 

B  =  arc  tan  (F  /F  ) 
y  x 

B.  THRUST  STAND  CALIBRATION 

The  thrust  stand  calibration  was  accomplished  according  to  procedures 
outlined  in  the  Gilmore  manual.  The  thrust  stand  was  calibrated  to  corres¬ 
pond  with  a  specific  set  of  test  conditions  which  were  a  side  thrust  force 
less  than  3,000  pounds  and  a  maximum  axial  thrust  of  30,000  pounds. 

These  values  were  derived  from  a  five-degree  thrust  vector  orientation 
and  an  axial  thrust  of  about  25,  000  pounds. 

The  calibration  was  performed  twice,  once  prior  to  propellant 
loading  and  verification  testing  and  once  after  the  firing.  The  first  cali¬ 
bration  was  performed  after  a  partial  trial  run  to  adjust  the  ranges  of  the 
digital  acquisition  system.  For  maximum  sensitivity,  thrust  stand  adjust¬ 
ments  were  made  before  the  calibration  to  insure  accurate  results.  The 
adjustments  were  very  minor,  being  a  simple  tightening  of  all  threaded 
connectors  in  the  side  force  load  trains.  The  first  calibration  was  done  to 
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verify  propel  functioning  and  provide  input  for  an  error  analysis.  The 
second  calibration  was  performed  to  verify  the  repeatability  of  the  first 
calibration.  The  calibration  equipment  (Figure  2)  consisted  of  a  highly 
accurate  load  cell  and  electrical  readout  box  installed  in  series  with  a 
two-way  hydraulic  ram.  A  single -calibration  ram  assembly  was  installed 
in  sequence  at  locations  Yjc»  Y Xjc>  and  Z£.  The  maximum 

loads  for  each  position  are  shown  in  Table  I. 


TABLE  I.  CALIBRATION  LOADS 


STATION 

MAXIMUM  CALIBRATION 
LOAD  (POUNDS) 

Xlc 

5000 

X2c 

5000 

Ylc 

5000 

Y2c 

5000 

Z 

40,  000 

c 

Side  thrust  calibration  loads  were  applied  in  20%  steps  from  zero  to 
maximum  tension  value,  through  zero  to  an  equal  compression  value  and 
then  returned  to  zero.  At  each  step  (registered  on  readout  box  and  con¬ 
trolled  by  the  hand-pump  hydraulic  power  supply)  the  force  measured  by 
each  load  cell  was  simultaneously  recorded  by  the  digital  data  acquisition 
system  The  "up  and  down"  run  was  repeated  three  times  at  each  cali¬ 
bration  location  before  progression  to  the  next. 

Data  obtained  from  the  first  calibration  attempt  led  to  the  following 
observations: 

1.  Z  axis  -  linear  and  reproducible,  no  hysteresis 
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X 

PLANE  CALIBRATION 

1 

Y 

PLANE  CALIBRATION 

Calib 

Force 

X1 

X2 

X3 

gfgJI 

Y1 

Y2 

X.c 

•°'5Xlc 

0 

5Xlc 

Ylc 

-0.  9592  Y. 

lc 

.0.  0408  X. 

1  c 

X2c 

-0  0204 

X, 

2c 

-0.  9592 

x, 

2c 

-0.  0204 

X, 

2c 

Y2c 

-0  0408  Y-v 

2c 

-0.  9592  YZc 

Figure  1.  Thrust  Stand  Dimensions  and 
Theoretical  Output  Force  Distribution 
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2.  X  axia  and  Y  axis  -  (a)  each  run  linear,  (b)  no  hyatereaia,  (c)  third 
runs  reproducible  when  compared  with  second  runs,  and  (d)  some  zero- 
shift  encountered  when  going  from  tension  to  compression  as  typified  by 

Fi«ure  4l  LOAD  (LBS) 

TENSION 


t  COUNTS 


COMPRESSION 

Figure  4.  Calibration  Counts  versus  Load 

Data  obtained  from  the  second  calibration  lead  to  the  following  ob¬ 
servations: 

1.  Z  axis  -  linear,  no  hysteresis,  reproducible  and  similar  to  the 
first  calibration. 

2.  Y  axis  -  linear,  no  hysteresis,  reproducible  and  similar  to  the 
first  calibration  except  the  compression  run.  This  run  was  nonlinear 
from  0  to  2000  pounds.  The  linearity  was  constant  from  2000  io  5000  lbs. 
The  Y^  axis  contributed  little  to  the  measurement  of  side  force  because 

the  side  force  was  generated  at  the  top  of  the  motor  near  the  nozzle.  From 
the  table  on  Figure  1,  it  can  be  seen  that  if  a  force  was  applied  at  Yjc, 
less  than  4.  2%  of  the  total  load  was  measured  by  the  Y^  load  cell.  During 
a  firing,  the  load  cell  would  measure  less  than  20%  of  the  total  load 
because  of  its  position.  The  nonlinearity  was  less  than  5%  from  the  first 
calibration  curve  from  0  to  1500  lbs  and  therefore,  the  nonlinearity  of 
Y^  contributed  about  1%  to  the  total  error.  Figure  5  shows  where  was 
nonlinear . 
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Figure  5.  Comparison  of  First  and  Second  Calibration,  Y  Load  Cell 
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3.  X  axis  -  each  run  linear,  no  hysteresis,  and  reproducible 
similar  to  the  first  calibration  except  the  axis  in  compression  was 
nonlinear  from  3000  to  5000  pounds.  The  maximum  nonlinearity  was  about 
10%  at  5000  pounds.  The  same  situation  exists  in  X,  as  in  except  there 
were  two  load  cells  on  top  instead  of  one  as  in  the  Y  direction.  measures 
less  than  15%  of  the  total  side  force,  therefore  the  error  contributed  to 
the  nonlinearity  of  X^  was  less  than  1,5%  at  5000  lbs  side  force.  Since  the 
maximum  anticipated  side  force  was  less  than  3000  lbs,  the  nonlinearity  will 
not  affect  side  force  measurement  in  this  range.  Figure  6  shows  the 
nonlinearity  of  X2. 

The  probable  reason  that  the  Xg  and  Y^  steps  in  compression  were 
nonlinear  at  certain  loads  was  because  the  load  cells  were  not  bolted  to  the 
structure  as  tight  as  possible.  About  5  months  elapsed  between  the  2 
calibrations.  During  this  time,  the  motor  was  subjected  to  the  following: 
a  fire  from  a  nearby  motor  malfunction,  significant  rain  and  wind,  as 
well  as  unknown  disturbances  by  normal  test  crew  operations  around  the 
motor . 

C.  THRUST  STAND  VERIFICATION 

The  purpose  of  the  verification  test  was  to  determine  how  accurately 
the  thrust  stand  functioned  as  a  measurement  system.  Each  of  the  load 
cells  had  been  calibrated  individually  in  the  laboratory  and  their  individual 
accuracies  determined.  The  effect  of  the  thrust  stand  geometry  on  mis¬ 
alignments  of  load  cells  and  indications  of  force  was  unknown,  so  the 
system  accuracy  was  determined  with  the  verification  test. 

Instead  of  conducting  a  firing,  which  normally  follows  a  stand  calibra¬ 
tion,  a  precisely  applied  load  simulated  the  thrust  and  side  force  (  in  the 
three  directions).  The  verification  loads  were  applied  near  the  predicted 
thrust  vector  action  point  on  the  Z  axis.  This  was  accomplished  by  the  use 
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of  a  highly  accurate  (secondary  standard)  load  cell  and  hydraulic  ram 
positioned  on  top  of  the  motor  (Figure  3).  Since  a  vertical  force  of 
40,  000  lbs  and  a  horizontal  force  of  5000  lbs  was  to  be  applied,  there  had 
to  be  a  structure  to  transmit  the  load  to  the  motor.  A  verification  struc¬ 
ture  (Figure  3)  was  fabricated  to  accomplish  this.  In  this  manner  known 
forces  were  applied  in  the  X,  Y  and  Z  directions. 

For  the  evaluation  of  the  X,  Y  and  Z  axes,  simultaneous  readings  of 
all  load  cells  were  recorded,  reduced  with  the  calibration  equations  and 
compared  to  the  known  values  of  loads  applied.  The  following  standard 
loads  were  applied  to  the  stand  in  20%  steps: 


X  direction 

5000  lbs. 

Y  direction 

5000  lbs. 

Z  direction 

40000  lbs. 

Since  the  verification  structure  was  not  completely  rigid,  it  deflected 
about  0.25  inches  when  a  vertical  force  of  40,000  lbs  was  applied.  When 
the  horizontal  force  of  5,  00^  lbs  was  applied  the  structure  deflected  in 
bending.  Due  to  this,  it  was  difficult  to  apply  standard  loads  exactly  in 
the  desired  plane.  If  the  load  was  misaligned  by  a  given  angle,  the  error 
in  indicated  applied  force  would  vary  as  the  cosine  of  the  angle.  It  would 
take  an  angle  of  8  degrees  to  produce  a  1%  error.  The  structure  was 
very  difficult  to  analyze  because  it  moved  in  a  complex  manner.  Thus, 
the  amount  of  error  that  was  attributed  to  the  movement  of  the  stand  was 
difficult  to  determine  exactly.  The  amount  of  misalignment  was  believed 
to  be  small  (less  than  1/2  degree)  because  the  surfaces  were  accurately 
machined  and  initially  aligned  with  accurate  machinist's  levels.  The 
evaluation  of  the  verification  results  confirmed  that  the  thrust  stand  pro¬ 
duced  acceptable  results  with  an  error  of  less  than  2%. 
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--DATA  THOM  1ST  CALIBRATION 

—  UATA  FROM  2ND  CALIBRATION 

THE  FIRST  SUBSCRIPT  REFERS 
TO  THE  LOAD  CELL  LOCATION 
AND  THE  SECOND  REFERS  TO 
THE  RUN  NUMBER 


•»•*•  »*•«»  •<»•«  -*«•«  -100J 


1  I  I  I  I  I  I  I 

BOOB  BB««  BOBO  IBBO  BOOB  ,000  10  0  0 

COUNTS 


rx“ 

/f^NON- 
rf  f  LINEAR 
*  PORTION 


POUNDS 


Figu  re  6.  Comparison  of  First  and  Second  Calibration,  X  Load  Cells 


SECTION  III  1 

ERROR  ANALYSIS 

A.  ANALYSIS 

The  following  errors  were  considered  In  this  analysis:  (1)  bias  error; 

(2)  errors  in  each  axis;  (3)  errors  in  the  side  force  calculation.  Figure 

7  shows  some  parameters  that  were  used  to  describe  the  forces  produced 

by  a  TVC  system,  but  they  are  not  discussed  in  regard  to  error.  The  i 

data  for  the  error  analysis  was  obtained  from  the  accuracy  verification.  r 

'( 

a 

i 

•f 

Z  ? 


Figure  7.  Vector  Description  of  Nozzle  Forces 

Each  error  includes  additively  several  preceding  errors.  For  example, 
the  error  in  the  side  force  contained  the  errors  from  the  X  and  Y  axis. 
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1 .  Bias  Error. 

The  data  from  the  verification  of  the  Gilmore  thrust  stand  shoved 
that  all  values  of  force  derived  from  the  calibration  and  data  reduction 
process  were  larger  than  applied  loads.  This  indicated  that  a  bias  error 
could  be  present.  This  bias  error  can  be  reduced  by  the  use  of  a  correction 
factor.  The  correction  factor  was  obtained  by  the  use  of  the  following 
formula: 


Specific  correction  factor  SCF  = 


input 


average  output 


where  input  =  specific  load  put  in  during  the  verification  test 

.  .  £  outputs  for  each  individual  input 

Ave-age  output  =  —  ■  - — - g - e — =  fJc 


X. 


There  were  6  different  outputs  obtained  from  each  input.  The  inputs 
for  the  X  and  Y  axes  were:  1000  lbs,  2000  lbs,  3000  ibs,  4000  lbs, 

5000  lbs.  The  inputs  for  the  Z  axes  were:  8000  lbs,  16,000  lbs,  24,000 
lbs,  32,  000  lbs,  4G, 000  lbs. 


The  total  correction  factor  (TCF)  was  obtained  by  averaging  the 
specific  correction  factors  for  each  axes.  The  total  correction  factors  for 
the  three  axes  are: 


X  axes  TCF  =  0.  9955 
Y  axes  TCF  =  0.  9876 
Z  axes  TCF  =  1.  0063 

More  specific  data  and  calculations  can  be  seen  in  the  Appendix,  Table  III. 
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2.  Errors  in  Each  Axis. 

The  error  in  each  axis  is  obtained  first  by  determining  the  standard 
deviation  c 


<y  = 


where  Xq  is  the  different  outputs  obtained  from  the  specific  input  and  //q 
is  the  average  of  all  the  outputs  as  defined  in  Part  A  (1)  of  the  error 
analysis,  (n  is  6:  n-1  is  used  because  of  the  small  sample  size  of  6). 


The  error  is  expressed  as  a  3c  error.  The  error  is  obtained  by 
dividing  3<r  by  the  average  outputs  for  each  axis. 

m  ♦  3  <r  dX 

%  point  error  =  tt-  =  — rr 

A  .X 


The  %  point  errors  are: 


X  axis  error  =  0.  44 
Y  axis  error  =  0.518% 

Z  axis  error  =  0.  18% 

More  specific  data  can  be  seen  in  Table  III  through  V  in  the  Appendix. 
3.  Errors  in  the  Magnitude  of  the  Calculated  Side  Force. 


The  side  force  is  given  by  the  relationship  Fg=^Fx^  +  and  the 


error  in  the  side  force  (dFg)  is  given  by  dFg  =  d(F^2 +  F  Let 


F  =  X;  F  =  Y;  dF  =  (2Xdx  +  2Ydy) 
3C  y  s 


■  T  *2  I  dr, ,  + 

(X2  +  Y2)  l^2'  I  (X2  +  Y2)  ^2 
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where 


X  =  X  side  force 
Y  =  Y  side  force 
dx  =  3<r  error  for  X  side  force 
dy  =  3<r  error  for  Y  side  force 


The  %  error 


XdX  +  Ydy _ i _  *  Xdx  +  YdY 

(X2  +  YV/2  (X2  +  Y2)1/2  X2  +  Y2 


The  error  in  the  side  force  can  be  made  up  of  many  different  combinations 
of  X  and  Y  loads  and  errors.  This  data  is  summarized  in  Table  VI  in  the 
appendix.  The  side  force  error  ranged  from  1.  14%  for  a  side  force  of 
1414  lbs  and  .  33%  for  7100  lbs.  In  general,  as  the  side  force  went  up,  the 
error  went  down. 


B.  DISCUSSION  OF  ERRORS 

The  error  analysis  for  the  Gilmore  thrust  stand  must  not  be  regarded 
as  the  ultimate  capability  of  the  Gilmore  stand.  This  error  analysis  was 
based  on  only  one  calibration  and  verification  before  the  firing  and  one 
calibration  after  the  firing.  The  calibration  after  the  firing  was  performed 
to  determine  if  the  stand  would  give  reproducible  results.  The  stand  did 
give  reproducible  results  as  mentioned  earlier  except  for  two  cases. 

Because  of  the  two  cases  which  did  not  reproduce  and  also  of  the  small 
sample  size  of  two,  the  confidence  limits  were  low.  The  test  is  very  expen¬ 
sive  (based  on  man-hours),  and  repeated  testing  to  increase  the  confidence 
limits  cannot  be  justified.  Continued  calibrations  and  analysis  during  the 
routine  use  of  the  thrust  stand  will  help  establish  repeatibility  of  the  cali¬ 
brations  and  confidence  limits  for  the  estimated  errors  of  the  thrust  stand. 


16 


C.  ERRORS  NEGLECTED 


Several  sources  of  errors  were  neglected.  These  include  angular 
misalignment  (cross  axis  effects)  and  dynamic  errors. 

1.  Angular  Misalignment. 


There  are  several  types  of  angular  misalignment  that  are  discussed 


below . 


a.  Motor  and  spider  are  misaligned  in  the  thrust  stand  while 
vertical  load  is  parallel  to  motor  centerline  (Figure  8). 


Figure  8.  Motor  and  Spider  Misalignment 

In  this  case  the  oflset  c.g.  creates  a  moment  about  0.  The  misaligned 
motor  creates  erroneous  readings  in  all  the  load  cells. 

b.  Motor  is  misaligned  in  thrust  stand  with  Z  spider  aligned 

correctly  while  vertical  load  is  parallel  to  centerline  (Figure  9). 

LOAD 


Figure  9.  Motor  Misalignment 
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In  this  case  the  load  does  not  act  through  0,  thus  creating  a  moment.  Also 
the  c.g.  creates  a  moment  about  0.  All  the  load  cell  readings  are  affected. 

c.  Motor  and  spider  arc  misaligned  in  thrust  stand  and  load  is 
misaligned  with  motor  (Figure  10). 


Figure  10.  Motor,  Spider  and  Load  Misalignment 

In  this  case  both  the  load  and  c.g.  create  a  moment  about  0.  All  the  load 
cell  readings  are  affected. 

d.  Only  the  load  is  misaligned  while  the  motor  is  assumed  to  be 
aligned  (Figure  11).  LOAD 


Figure  11.  Load  Misalignment 
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The  misaligned  load  would  create  errors  in  all  load  cell  readings.  The 
error  in  the  Z  direction  equals  (1*  cos  (1)  and  the  error  in  the  X  and  Y  direc¬ 
tion  equals  s  in  p  (Figure  12).  Only  the  last  type  of  misalignment  is  dis¬ 
cussed,  because  in  the  other  three  cases  it  was  difficult  to  determine  the 
amount  of  misalignment  or  because  there  were  combinations  of  two  types 
of  misalignment. 


Figure  12.  Side  Load 


A  known  vertical  load  applied  to  the  motor  at  an  8 -degree  misalignment 
would  cause  an  error  in  the  Z  axis  of  1%-  In  contrast,  an  angular  misalign¬ 
ment  cf  only  35  minutes  of  arc  in  the  vertical  load  would  cause  the  X  and  Y 
load  cells  to  read  out  1%  of  the  vertical  load.  This  is  clarified  by  Table  II. 


TABLE  II.  CROSS -AXIS  EFFECTS 


e 

100  Sinp  =  % 

Cos  p 

100  (1  -cosp)  =  % 

1  min 

0. 00029 

1. 00000 

0. 00000 

10  min 

0. 00290 

1. 00000 

0. 00000 

35  min 

0. 01018 

0.99995 

0. 00005 

1  degree 

9  min 

0. 02007 

0. 99980 

0. 00020 

2  degrees 

52  min 

0. 05001 

0. 99876 

0. 00124 

8  degrees 

6  min 

0.  14090 

0. 99002 

0. 00998 
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A  cross -axis  effect  is  the  reading  of  a  load  cell  which  is  not  exactly 
perpendicular  to  the  applied  load  action  line.  Even  with  a  small  angular 
misalignment,  some  cross -axis  effects  exist. 

2.  Dynamic  Errors. 

The  response  (in  terms  of  deflection)  of  a  linear  spring-mass 
damper  vibrational  system,  when  excited  by  a  time -varying  force,  exhibits 
both  time  and  magnitude  "errors"  when  the  dynamic  deflection  is  compared 
to  the  static  deflection  under  the  same  load.  In  the  case  of  a  rocket  motor 
vertical  thrust  stand,  additional  complexity  results  from  the  varying  mass 
of  the  motor  as  propellant  is  consumed.  Since  load  cell  output  is  a  direct 
function  of  load  cell  deflection,  "errors"  in  deflection  result  directly  in 
errors  in  the  load  cell  output.  While  it  is  theoretically  possible  to  compute 
the  input  (applied  load)  from  the  dynamic  stand  response  (load  cell  reading) 
the  multiple  degrees  of  freedom  present,  as  well  as  the  nonlinear  spring 
constants  of  *v  vibrational  elements,  make  this  approach  a  practical 
impo,.„ibilit, . 

The  practical  alternative  is  to  avoid  exciting  the  system  at  frequencies 
near  its  nature  frequencies ,  particularly  the  lowest  natural  frequency.  In 
terms  of  test  procedures,  this  requires  that  axial  and  side  thrust  be  varied 
at  fairly  low  cyclic  rates  (on  the  order  of  2  to  5  cycles  per  second).  The 
allowable  cyclic  rates  could  be  increased  by  "stiffening  up"  the  stand  (i.e. 
by  eliminating  flexures  and  using  higher -capacity  load  cells),  but  this  in 
turn  would  result  in  poorer  static -load  measurement  accuracy. 
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SECTION  IV 


RESULTS  AND  CONCLUSIONS 

The  error  analysis  performed  on  the  Gilmore  50,  000-lb  thrust  stand 
calibration  data  established  that  there  was  a  .  58%,  1. 24%,  and  >.625%  of 
point  bias  in  the  X,  Y,  and  Z  axes  respectively.  This  was  partially 
accounted  for  by  the  use  of  a  correction  factor.  The  aide  force  magnitude 
error  ranged  from  about  1.  1%  to  .  3%  of  point  throughout  the  measurement 
range.  This  did  not  take  into  account  cross -axis  effects  introduced  by  a 
single  vertical  load. 

The  major  problem  was  adequately  defining  the  errors  in  measured 
forces  due  to  thrust  stand  interactions  (cross -axis  effects),  changing  motor 
weight  and  axial  center  of  gravity  location  as  a  function  of  time.  The  com¬ 
puter  program  that  is  described  in  the  appendix  will  calculate  the  decrease 
in  motor  weight  with  time,  thus  reducing  the  error  associated  with  this 
factor.  The  error  in  the  aide  force  was  small,  less  than  1.  1%  of  point, 
considering  the  inaccuracy  of  the  stand  alignment.  If  the  angular  misalign¬ 
ment  of  the  stand  was  known,  an  error  could  be  calculated  to  take  misalign¬ 
ment  into  account.  Since  it  was  not  known,  the  only  errors  not  taken  into 
account  were  dynamic  and  misalignment  errors. 


SECTION  V 


RECOMMENDATIONS 

The  Gilmore  six-component  thrust  stand  gave  acceptable  results. 

Since  dynamic  errors  were  not  taken  into  account,  the  error  associated 
with  this  source  was  unknown.  Also  the  amount  of  angular  misalignment 
was  unknown.  Thus,  if  a  more  accurate  estimate  of  the  error  is  needed, 
a  dynamic  error  analysis  and  an  accurate  stand  alignment  (possibly  optical) 
should  be  performed.  Also  if  this  error  is  to  be  accepted  with  high  con¬ 
fidence,  at  least  four  more  calibrations  and  verifications  should  be  per¬ 
formed  on  the  test  stand. 


22 


I 


$ 


APPENDIX 


CALCULATIONS 
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APPENDIX 


CALCULATIONS 


A.  ERROR  ANALYSIS 

The  following  errors  were  considered:  (1)  error*  in  each  exit; 

(2)  error  in  the  aide  force. 

1.  Error  in  Each  Axi». 

The  data  from  the  calibration  of  the  Gilmore  thrust  stand  showed 
that  all  values  of  force  applied  were  larger  than  the  applied  loads.  Thus, 
there  was  a  bias  error.  This  bias  error  can  be  reduced  by  using  a  correc¬ 
tion  factor. 

Sample  Calculation: 


OUTPUT  (lbs) 


X  INPUT 

RUN  1 

RUN  2 

RUN  3 

MODE 

1000 

1016.79 

1005. 75 

1003. 43 

Increasing  Load 

2000 

5000 

io:oo 

1005.60 

1010.  19 

1003.  08 

Decreasing  Load 

Tae  average  for  all  the  1000-lb  inputs  was  1007.47. 


Correction  factor  c.f.  =  1000.00/1007.47  =  .9925 
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In  the  same  way  correction  factors  were  obtained  or  other  inputs.  These 
are  shown  in  Tables  III  through  V. 

The  correction  factor  for  the  X  axis  is:  0.9955 
The  correction  factor  tor  the  Y  axis  is:  0.9876 
The  correction  tactor  for  the  Z  axis  is:  1.0063 

These  were  obtained  by  averaging  the  correction  factors  for  each  axis. 

The  3<r  error  is  obtained  by  first  averaging  values  for  each  input  X, 

T.  Z.  From  this,  o-  is  obtained  by  the  formula  <r  =  /lifX^-X)^  where  N  is 

V 

the  number  of  indicated  values.  N-l  is  used  because  of  the  small  sample 
size  of  6. 

Sample  Calculations: 

x  =  1007.47 

x,  -  x  -  1005.60  -  1007.47  =  -  1.  87  lbs 
x*  .  x  =  1005. 75  -  1007.47  t  -  1. 72  lbs 


x^  -  x  =  1003.  08  -  1007.  47  =  -4.  39  lbs 
=  (1.87)2  +  (1.72)2  +  .  .  +  (4. 39)2 
or  1000  =  3.  515:%  point  error  -  3(r/Xi  -  10*  55/100C  =  1  •  05% 
<r  4000  -  4.51:%  point  error  =  3<r/Xi  =  13.53/*i000  =  .  34% 

The  rest  of  the  values  are  shown  in  Tables  III  through  V. 
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TABLE  III.  LOAD  CELL  EVALUATION  FOR  X 
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uracy  of  Primary  Standard  -  99.  99% 
uracy  of  Secondary  Standard  -  99.95% 


TABLE  IV.  LOAD  CELL  EVALUATION  FOR  Y 


TABLE  V.  LOAD  CELL  EVALUATION  FOR  Z 


Accuracy  of  Primary  Standard  =  99.  99% 
Accuracy  of  Secondary  Standard  =  99.  95% 


Error  in  the  Side  Force. 


The  side  force  F 


2  2 
+  F 

y 


The  error  in  the  side  force  =  dF 


-  d(F*2  *  Fv2) 


2\l/2 


Let  F  = 
x 

F  = 


X 

Y 


dF,  .  <2Xd*  4  ZYdyj-ip-i-jpj 


dF 


Xdx  +  Ydy 

-  i\tn 


(x2  4  Y2) 


Sample  calculation 


dFg  for  X  =  1000  lbs,  Y  =  1000  lbs 


3a 1000-x  =  10. 55  lbs 


3a  1000-y  =  12.  30  lbs 


The  above  values  are  obtained  from  Tables  III  through  V. 


dF 


7b  error 


(x2  4  Y 2)'/2 


Xdx  +  Ydy  1 

(x2  +  Y2)1'2  '  {?  4  Y2),/2 


Xdx  +  Ydy 
— — 2 
X  +  Y 
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For  the  above  values 


1000(10. 551  4  1000(12.  30) 

10002  4  10002 


10550  +  12300 

2, 000, 000 


22850 

2x1? 


1.  14% 


The  rest  of  the  values  are  shown  in  Table  VI  and  Figure  13. 
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TABLE  VI.  SIDE  FORCE  BIAS  AND  ERROR 


X  LOAD 
X 

Y  LOAD 
Y 

SIDE  FORCE 
Fs 

SIDE  FORCE  ERROR 

dF  /F  % 

8  8 

SIDE  FORCE  BIAS 
% 

1000 

1000 

1414 

1.  14 

.  82 

1000 

2000 

2236 

.  55 

.  74 

1000 

3000 

3163 

.  M 

.  1? 

1000 

4000 

4123 

.  22 

.  14 

1000 

5000 

5099 

.  33 

.  19 

2000 

1000 

22  36 

.  85 

.82 

2000 

2000 

2828 

.  59 

.  77 

2000 

3000 

3605 

.  54 

1. 08 

2000 

4000 

4472 

.29 

1.  33 

2000 

5000 

5  385 

.  37 

1.77 

3000 

1000 

3163 

.49 

.71 

3000 

2000 

3605 

.42 

.  71 

3000 

3000 

4243 

.  43 

.95 

3000 

4000 

5000 

.  26 

1.  19 

3000 

5000 

5831 

.  33 

1.  58 

4000 

1000 

4123 

.  39 

.  33 

4000 

2000 

4472 

.  35 

.  39 

4000 

3000 

5000 

.  38 

.62 

4000 

4000 

5657 

.  25 

.88 

4000 

5000 

6403 

.  31 

1.29 

5000 

1000 

5099 

.  39 

.  29 

5000 

2000 

5  385 

.  37 

.  19 

5000 

3000 

5831 

.  38 

.  07 

5000 

4000 

6403 

.  28 

.  36 

5000 

5000 

7071 

.  33 

.  79 

All  of  the  values  of  X  and  Y  Load  were  high  except  for  one.  Thus,  all 
values  of  side  force  would  be  high.  Thus  the  correction  factor  for  the  side 
force  would  be  less  than  1,0,  i.e.,  Correction  Factor  =  1.0  -  Bias  Error. 
Figure  9  shows,  in  general,  as  the  magnitude  of  the  side  force  increases, 
the  error  in  the  magnitude  of  the  side  force  decreases  as  shown  by  the 
points.  The  amount  of  bias  is  a  function  of  the  X  and  Y  Loads  for  the  same 
side  force.  Thus,  a  side  force  of  3605  lbs  can  have  a  bias  ranging  from 
1.42%  to  0.  71%. 
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8000 
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Figure  13.  Side  Force 


Summing  forces  and  talcing  moments: 

' i  T  "2  T  '*3 


ZF„  =  0F.  =  X,+X2+X3  1)  i:Mx  =  FyZ  -  FZY  -  Y|  D  -  0 


SFy  =  0  Fy  =  Y1  +  Y2 


SF  =  0  F  =  Z 
z  * 


2)  2My  =  FzX  -  FXZ  +  (Xj  +  X3)  D  -  0 


3)  ZMx  =  FyX  -  FxY  +  (Xj  +  Y3)  W  =  0 


From  1) 


^  F2Y  +  YjD 
Z  =  - P - 
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Sub  1) into  2 


-F  F  Y  -  F  Y.T5 

-  T - 2-L_  +  FtX  +  (XA  +  X3)T5  s  0 


or  4) 


_  F  F  Y  +  F  Y,D  w 

F  X  *  ■*-*  V  *  *  -  (X1  +  x3)T5  s  0 

*  y 


or  4) 


F„F_Y  +  F„Y,D  -  (XL  +  X3>15F 


z  T  x‘  1~ 

F  F 

y  * 


4)  into  3 


tFzY  4  TxYxD  -  (Xj  +  X3)DFy  -  F^Y  4  (Xt  -  X3)W 


F  Y,D  -  (X,  4  X,)  DF  +  (Y.  -  Y ,)  W 
s  x  1 _  1  3'  y  1 _ 3  q 

F 

z 

The  two  terms  with  unknowns  will  cancel,  thus  there  is  no  solution  for  Y. 


Try  the  same  approach  but  with  different  values 


.  I 


From  2) 


F  X  =  F  2  -  (X,  +  X,)D 
z  x  1  3 


FxZ  -  (Xj  +  X3)D 


"1 


into  3 


F  F  Z 

-  y  * 


F  (Xj  ♦  X3)D 


FXY  -t  (Xx  -  X3)W  *  0 


or 


into  1) 


_  ryrx7  -  Fy(X,  +  x3)d  .  (x,  -  x^f, 


rxFyZ  •  F,F,2  •  V*1  >  X,)D  ,  (X,  -  X3)WFy  ^ 

F* 

Again  the  two  terms  with  unknowns  cancel  each  other.  Thus,  the  values 
for  two  of  the  components  cannot  be  computed.  The  only  one  that  can  be 
determined  is  Z. 

Thus,  taking  out  the  components  X  and  Y,  Z  can  be  computed. 

EM  =  Y,D  -  F  Z  *  0 
x  l  y 

EM  »  (X,  4  X,(D  -  F  Z  =  0 
y  1  3'  •  x 

EM  =  (X,  +  X,)W  =  0 

S  1  3' 

Solving  for  Z 

Z  =  Yl°  or  z  ,  (X>tX3)15 

y,.y2  (x1+x2  +  x3) 

These  stand  equations  give  the  location  of  point  of  application  of  the  thrust 
vector. 
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C.  COMPUTER  PROGRAM 

A  computer  program  was  written  to  provide  the  following  item*: 

1.  Magnitude  and  direction  of  side  force  vector. 

2.  Axial  thrust  corrected  to  given  mass  flow  rate. 

3.  Ratio  of  resolved  aide  force  to  axial  thrust. 

4.  Amplification  factor. 

During  a  firing,  propellant  is  being  consumed,  and  the  center  of  gravity 
is  constantly  changing.  The  weight  of  the  expended  propellant  exhausted  was 
determined  by  a  computer  program.  The  correct  value  of  thrust  was 
obtained  in  the  following  manner: 

Weight  of  propellant  exhausted  at  anytime  t  =  Wpt  (lbs)  -  burn  rate 
lbs/sec  X  incremented  time  T  (sec)  =  (lbs) 

Force  at  anytime  read  by  load  cells:  t  =  F{ 

Total  weight  of  propellant  at  start  =  Fp 

Correct  Z  force  =  F  z  -  (Fp  -  Wpt)  =  corrected  thrust  at  any  time  t. 

The  thrust  values  obtained  from  the  load  cells  are  corrected  and  the  true 
thrust  is  plotted  on  graph  paper  as  thrust  versus  time. 
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~*Tnis  report  describes  the  calibration  and  an  error  analysis  of  the  6-component 
50,  000-pound  Gilmore  Rocket  Motor  Thrust  Stand  located  on  Pad  2,  Air  Force 
Rocket  Propulsion  Laboratory,  Solid  Test  Area  (1-32;.  The  purpose  of  this 
report  is  to  establish  stand  accuracy  and  reproducibility  for  a  calibrated  static 
test  loading. 


Errors  considered  primarily  resulted  from  thrust  stand  misalignment, 
propellart  mass  flow  cf  the  motor,  and  cross-axis  interactions.  Dynamic  effects 
suiting  from  vibration  were  not  included  in  the  detailed  analysis  but  are  briefly 
icussed. 
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